Nuclear Matter and Neutron stars in a Parity Doublet Model 
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We investigate the properties of isospin-symmetric nuclear matter and neutron stars in a chiral 
model approach adopting the SU(2) parity doublet formulation. This ansatz explicitly incorporates 
chiral symmetry restoration with the limit of degenerate masses of the nucleons and their parity 
partners. Instead of searching for an optimized parameter set we explore the general parameter 
dependence of nuclear matter and star properties in the model. We are able to get a good description 
of ground state nuclear matter as well as large values of mass for neutron stars in agreement with 
observation. 



INTRODUCTION 

The use of effective hadronic models based on chiral 
symmetry has been quite successful in describing nuclear 
matter, finite nuclei and neutron stars (see e.g. 
While most of the approaches are based on variations of 
the linear sigma model or its non-linear realization, work 
on an alternative ansatz starting from the parity-doublet 
model (0,0,1) has been quite limited. In this approach, 
the scalar sigma meson serves to split the nucleon and 
its parity partner, while in the chirally restored phase 
both baryonic states become degenerate but not mass- 
less. This situation is quite analogous to sigma models 
with hidden gauge symmetry [B[ where the masses of vec- 
tor and axial vector mesons are split in a similar way to 
attain degenerate (but non-zero) masses in the chirally 
restored phase similar to the pion and sigma mesons. It 
has only quite recently been shown [f| that with this 
approach one can obtain saturated nuclear matter with 
relatively high nuclear compressibilities with a simple po- 
tential for the scalar field and a cross-over transition to 
chiral restoration at large densities. In our investigation 
we calculate symmetric nuclear matter and neutron star 
properties within this framework. In order to get a bet- 
ter understanding of the parameters involved and not just 
to fix an arbitrary fine-tuned parameter set we vary the 
parameters over a broad range of values and calculate 
nuclear matter properties and star masses. 

The parameter scan shows that for quite a wide range 
of values a satisfactory reproduction of nuclear matter 
binding energy and saturation density is possible. Fur- 
ther demanding that the nuclear matter compressibility 
lies within a range of reasonable values (assumed to be 
between 100 and 400 MeV in the following), restricts the 
parameter space significantly but still allows for a range 
of values that yield a quantitatively good description of 
nuclear matter ground state properties in contrast to the 
results in 0]. In order to study isospin asymmetric nu- 
clear matter the rho meson is included in the Lagrangian. 



Its coupling to the nucleons is fixed by the phenomenolog- 
ical value for the asymmetry energy of about 32 MeV. By 
further requiring that the maximum mass of the neutron 
star, determined by solving the Tolman-Oppenheimer- 
Volkov (TOV) equations, should have a mass around 1.9 
M s olar, we find additional restrictions for the allowed 
model parameters. 

The article is organized as follows. The parity doublet 
model is briefly introduced in the following section. We 
then apply the equations to symmetric nuclear matter 
and to isospin-asymmetric matter as it is present in neu- 
tron stars. The properties of the model are first tested 
to fit saturation quantities and then to reproduce ob- 
servational properties such as the maximum mass of the 
neutron star. Finally, we discuss the results of this in- 
vestigation and add some remarks on further work to be 
done. 



MODEL DESCRIPTION 

In the parity doublet model the nucleons and their par- 
ity partners belongto the same multiplet and so are true 
chiral partners Under SU L {2) x SU(2) R trans- 

formations L and R, the two nucleon fields ipi and tp2 
transform as: 
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This allows for a chirally invariant mass, mg: 
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The full chiral Lagrangian reads: 

C = + 02^02 

+ m (027501 - "017502) 

+ ai>x (cr + i-f 5 r • 7r) 0>i + fr0 2 (c - «75T • tt) 02 

- .9^0i7m lj ' 1 V'i - 9uji>2l^ l i> 2 

- 9p^\lp T ■ ^01 - 9p^2lpT ■ ^02 

+ £m , (4) 

where a, 6, g w and g p are the coupling constants of the 
mesons fields (cr, tt, uj and p) to the baryons 0i and 02. 
Note that we assume the same vector coupling strength 
for both parity partners, as it would result from min- 
imal coupling of the vector mesons. The mesonic La- 
grangian C m contains the kinetic terms of the different 
meson species, and potentials for the scalar and vector 
fields. The potential for the spin zero fields is the same 
as in the ordinary SU(2) linear sigma model. Kinetic and 
potential terms are added for an isoscalar vector meson, 
u>, as in the a-to model of nuclear matter [3| as well as 
for the isovector vector meson, p: 

Cm = + ^ M 7pd% - 
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where uj pv = d^ujv - d v uj p and p M „ = d^p v - d v p p rep- 



resent the field strength tensors of the vector fields. As 
usual, the parameters A, p and e can be related to the 
sigma and pion masses, and the pion decay constant, in 
vacuum: 

A = 



2 erg 



(6) 



with m^ = 138 MeV, /„. = 93 MeV and the vacuum 
expectation value of the sigma field cro = fir- Since the 
mass of the a meson in the vacuum can not be fixed pre- 
cisely by experiment, we will treat it as a free parameter. 
The vacuum mass of the u> field is m u = 783 MeV and 
the vacuum mass of the p field is m p = 761 MeV. 

To investigate the properties of dense nuclear matter 
and the chiral phase transition at zero temperature, we 
adopt the mean-field approximation Q. The fluctua- 
tions around constant vacuum expectation values of the 
mesonic field operators are neglected, while the nuclc- 
ons are treated as quantum-mechanical one-particle op- 
erators. Only the time-like component of the isoscalar 



vector meson uj = luq and the time-like third component 
of the isovector vector meson p = po 3 survive, assuming 
homogeneous and isotropic infinite nuclear matter. Ad- 
ditionally, parity conservation demands tv — 0. The mass 
eigenstates for the parity doubled nucleons, the N + and 
N- are determined by diagonalizing the mass matrix, 
Eq. ([3]), for 0i and 02- Writing the coupling constants a 
and b as a function of M/v + (939MeV), Mm_ and cro, the 
effective masses of the baryons are given by: 
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If chiral symmetry is completely restored, i.e. a = 
0, the two nucleonic parity states become degenerate in 
mass with = — niQ. Thus, for large values of 

mo, the sigma field mainly serves to split the nucleon and 
its parity partner nucleon masses and contributes less to 
the (dynamical) mass of the nucleon Q. Independent 
of the value for mo, in the vacuum, i.e. for a — cro, the 
baryon effective masses reproduce their measured values. 

The thermodynamic potential per unit volume is: 



— = -c 



M 



^ (2tt) 3 



d 3 k(E*(k)-p*) , (8) 



where i G {V + ,AV} denotes the nucleon type, ji is 



the fermionic degeneracy, E*{k) = ^k 2 +M* z the en- 
ergy, and p* = m - g^uio - g p p 03 T 3 /2 = \J k 2 F + M* 2 
the corresponding effective chemical potential. The sin- 
gle particle energy of each parity partner i is given by 
Ei(k) = E*(k) + g^ujQ + g p p Q3 T 3 /2. 

The mean meson fields a, uj and p are determined by 
extremizing the thermodynamic potential £l/V: 
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The scalar density p si and the baryon density pi for each 
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chiral partner are given by the usual expressions: 
" feF * d 3 k M* 
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The basic nuclear matter saturation properties we im- 
pose can be formulated in the following way. At the point 
of vanishing pressure we have 



E/A-M N = -16.0±0.5MeV , 
po = 0.15 ±0.015 fm~ 3 



(11) 



Note that in contrast to former work in [6| , here we allow 
for a reasonable range of values for the different observ- 
ablcs. 
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FIG. 1: Value of the nucleon vector coupling constant gjvu vs 
mass parameter mo. 



RESULTS 
Nuclear Matter 

First we apply the model to symmetric nuclear matter 
and determine the ranges of the parameters that repro- 
duce the saturation properties for binding energy per nu- 
cleon and baryon density as given by eq. 1111 In principle 
the free parameters of the model are the coupling con- 
stant of the vector meson to the nucleons gNu>, the mass 
of the chiral partners M^_ (as the identification with 
a specific experimentally known resonance is not com- 
pletely unambiguous), the bare mass too, the coupling 
constant of the fourth order self interactions of the vec- 
torial fields <?4 and the mass of the scalar meson m a as 
already described in [6|. However, to keep things simple, 
we will for most of this work restrict our analysis to pa- 
rameter sets with 54 = and Mm_ = 1535 MeV. Where 
necessary or interesting we will discuss the influence of 
finite values of 54 or smaller values of Mm_ , as was also 
done in The other parameters are varied over all 
possible values. In Fig. |T]) and Fig. ((2]) we show the 
values of the corresponding parameters vs mass parame- 
ter mo, for which the imposed nuclear matter saturation 
properties are fulfilled. In both figures, the light grey 
bands represent the allowed regions without any further 
restriction, while the dark grey bands depict the param- 
eter sets, which in addition yield a reasonable nuclear 
matter compressibility K = 9f|j- = 9^0 = ^Pi^T at 
saturation between 100 and 400 MeV. As one can see, 
successful fits without restricting the compressibility are 
possible for 300 < too < 950 MeV. This region is larger 
than that found in [6( because no fine-tuning to fixed 
value for the nuclear matter saturation properties has 



been done here. The corresponding nucleon-vector cou- 
pling constant g^^ as well as the value of the sigma mass 
in vacuum to ct decrease with increasing values of too, as 
has also been observed in [9]. 
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FIG. 2: Mass of sigma meson in vacuum vs mass parameter 
mo. 

Demanding the nuclear matter compressibility to be 
in "reasonable ranges" restricts the possible parameter 
values. For our choice of 100 MeV < K < 400 MeV, 
we see that only the fits with quite large too values (> 
800 MeV) survive. 

This can be seen in more detail from Fig. (J3J. For 
smaller values of too the nuclear matter compressibility 
increases very fast beyond any acceptable value. 



4 



5000 




300 400 500 600 700 800 900 
m [MeV] 



FIG. 3: Nuclear compressibility at saturation vs mass param- 
eter mo. 



order self interactions of the vectorial fields produce more 
massive stars, being in agreement with phenomenology. 
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FIG. 4: Mass of neutron stars vs mass parameter mo. 



Asymmetric Matter and Neutron Stars 

In order to study neutron star properties we calculate 
the equation of state (EOS) of charge neutral matter in- 
cluding leptons. Due to the different densities of neu- 
trons and protons, the expectation value of the vector- 
isovector meson p does not vanish. The corresponding 
coupling constant is determined such that an asymmetry 
energy at saturation of E sym — 32.5 MeV is obtained. 
Neutron star masses and radii are determined by solv- 
ing the Tolman-Oppenheimer-Volkov equations [H, [Io| . 
The maximum star mass for a given equation of state 
is shown as a function of mo in Fig. (T?|) for the cases 
that have a compressibility at saturation K between 100 
and 400 MeV. As can be seen, the maximum mass of 
the star decreases with increasing tuq. This results from 
the fact that for larger bare mass term mo, medium ef- 
fects are reduced which leads to a smaller compressibility 
as can be seen in fig. [3] and reduced repulsion at high 
densities and consequently lower star masses. Thus, de- 
manding the compressibility to be in a reasonable range, 
especially not to be too big and in addition, to demand 
to have a maximum star mass above a phenomenologi- 
cally acceptable value of at least 1.4 to 1.6 solar masses, 
restricts the value of mo to be somewhere between 900 
and 870 MeV. Thus the initially broad band of possible 
parameter sets is strongly narrowed by these additional 
observables. Note, that these values of mo are quite large, 
which means that the mass generation is very much dif- 
ferent from a standard linear a model, where the mass 
is more or less exclusively generated by the dynamical 
symmetry breaking, while now most of the mass is gen- 
erated by the parity-mixing mass term. The results for 
neutron stars show that keeping a high mass for the chi- 
ral partners and a small coupling constant of the fourth 



In figure OH the central densities of the neutron stars 
and the threshold densities above which the A_ states 
are populated, respectively, are shown. The maximum 
mass stars have a typical central energy density of around 
lfm~ 3 which for the cases at hand with 54 = and 
Mn_ = 1535 MeV is always smaller than the critical 
density to create chiral partners and in this case no JV_ 
are present in the star as can be seen in the particle den- 
sities shown in Fig. [6] However as has already been 
observed in Q , decreasing the mass of the iV_ , reduces 
the threshold density and thus can allow for the popula- 
tion of iV_ states in neutron stars. This can for example 
be seen in Fig. for M N _ = 1379 MeV, where the 
chiral partners become populated already slightly below 
Pb = 1 fm~ 3 . The mass-radius relation for the two cases 
is shown in Fig. (jSJ. Besides the EOS for the dense part 
of the star, the EOS for an outer crust, an inner crust 
and an atmosphere have been added 11 1 . 



As can be seen, the maximum mass of the star is higher 
for M N _ = 1535 MeV (M max = 1.87 M Q ) compared 
with the case for M N _ = 1379 MeV {M max = 1.77 M c ), 
as expected. However, both masses are around 1.9 M Q , 
value that is in agreement with the massive stars observed 
lately with the heaviest one so far having a mass of M = 
2.ll° 'J M 0. 
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FIG. 5: Central densities of neutron stars vs mass parameter 
mo and critical density for chiral restoration or density for 
appearance of N-, respectively. 
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FIG. 7: Population of the star vs baryon density for mo = 855 
MeV and M N _ = 1379 MeV. The central density of the 
respective star would be pcent = 1.04/m -3 . 
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FIG. 6: Population of the star vs baryon density for mo = 843 
MeV and Mn_ = 1535 MeV. The central density of the 
respective star would be p ce nt = 0.96/m -3 . 



CONCLUSIONS 

We applied the parity doublet model to nuclear and 
neutron star matter. We were able to achieve a quan- 
titatively good description of saturated nuclear mat- 
ter as well as to obtain maximum masses of neutron 
stars that are in agreement with recent observations. 
In our approach we explored a wide range of model 
parameters that reproduce saturation properties of nu- 
clear matter for a window of values for the density 
(po = 0.15/to -3 ± 10%) and binding energy per nucleon 
(E/A - M N = -16 ± 0.5MeV). With that we could 
study the dependence of results on the free parameters 
of the model: the coupling constant of the vector meson 
to the nucleons gjv^, the bare mass toq and the mass 




10 20 30 

R [Km] 



FIG. 8: Mass radius diagram for different chiral partners 
masses. 

of the scalar meson m a . Initially the mass of the chi- 
ral partners is assumed to be M^_ — 1535MeV, cor- 
responding to one of the likely candidates listed in the 
particle data book. No vector meson self-interactions 
have been taken into account as they generally lead to 
very soft equations of state and small masses of neutron 
stars that are in disagreement with observation. Besides 
that, if we take into account phenomenological values 
of compressibility (allowing for a somewhat wider range 
lOOMeV < K < 400MeV) we find that the mixing pa- 
rameter mo has to be tuq > 800MeV so that the two 
different couplings of the nucleon to the scalar field can- 
cel each other to a large degree. 

In a second step the asymmetry between protons and 
neutrons has been included demanding a charge neutral 
state (including the electron) and the p meson coupling 
constant is determined to reproduce the phenomenologi- 
cal value for the asymmetry energy at saturation. If we 
consider a high mass for the nucleon chiral partners Mn_ , 
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the critical density above which those particles appear in 
the star is also high, and in the case of Mm_ = 1535MeV 
they do not appear in the star. The value of the max- 
imum mass is M max = 1.87M . Demanding that the 
M max should be at least 1.4 solar masses restricts the 
value of the explicit mass parameter mo to be less than 
900 MeV. Decreasing the mass of the chiral partners by 
about 150MeV is enough to decrease the threshold den- 
sity of the chiral partners for them to be present in the 
star, which reduces the maximum mass of the star mod- 
erately, in the case of M N _ = 1379MeV to M max = 1.77 
M Q . An interesting extension of this work would be the 
generalization to a Flavor-SU(3) approach. Work in this 
direction is in progress. 
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